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Abstract. The complete 0(a) QED initial state (IS), final state (FS) and initial-final state (IFS) interfe- 
rence corrections to the process e + e~ — > tt + tt~ axe presented. Analytic formulae are given for the virtual 
and for the real photon corrections. The total cross section (a), the pion angular distribution (da jd cos 0) 
and the tt + -k~ invariant mass distribution (da/ds') are investigated in the regime of experimentally realistic 
kinematical cuts. It is shown that in addition to the full 0(a) corrections also the 0(a 2 ) and leading log 
0(o?) photonic corrections as well as the contributions from IS e + e~ pair production have to be taken into 
account if at least per cent accuracy is required. For the data analysis we focus on an inclusive treatment of 
all photons. The theoretical error concerning our treatment of radiative corrections is then estimated to be 
2 per mill for both the measurement of the total cross section and the tt + -k~ invariant mass distribution. 
In addition we discuss the model uncertainty due to the pion substructure. Altogether the precision of the 
theoretical prediction matches the requirements of low energy e + e~ experiments like the ones going on at 
DA<£NE or VEPP-2M. 



1 Introduction 



Tests of the Standard Model (SM) as well as establishing 
possible new physics deviations from it crucially depend 
on our ability to make precise predictions. This requires in 
the first place a precisely known set of independent input 
parameters, like the fine structure constant a, the Fermi 
constant G M and the Z boson mass Mz- In fact screen- 
ing by vacuum polarization (VP) leads to an energy-scale 
dependent running electromagnetic coupling constant 



The precise definition of Aa(s) reads 

Aa(s) = -47raRe [Ilfo) - 77^(0)] , 



(2) 



where n'(s) is the photon vacuum polarization function 



(3) 



a(s) 



1 - Aa(s) 



G) 



of which the precise knowledge is crucial for elcctroweak 
precision physics. This effective coupling is sensitive to 
vacuum polarization effects, with about equal contribu- 
tions from leptons and quarks, causing the shift Aa(s), 
which is the sum of the lepton (e, /i, r) contributions and 
the contribution from the 5 light quark flavors (u, d, s, c, b): 

Aa(s) = Aai C p(s) + Aa^ d (s). At higher energies y/s > 
Mz also the heavier charged particles, the W and the top 
quark contribute. 



^yd 4 x e ^(0|T J ^(x)J: m (0)|0) 

= -(? 2 .r-9YK(? 2 ) 

and j£ m (x) is the electromagnetic current. 



Lcptonic contributions can be calculated perturbatively. 
However, due to the non-perturbativc behavior of the strong 
interaction at low energies, perturbative QCD only allows 
us to calculate the high energy tail of the hadronic (quark) 
contributions. Thus the main difficulty to determine the 
relationship between the low energy fine structure con- 
stant and the effective one at higher energies is the accu- 
rate determination of the non-perturbative contributions 
from low energy hadronic vacuum polarization insertions 
into the photon propagator. 

A way which allows us to do this is the precise mea- 
surement of low energy hadronic cross sections <7had(s) = 
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Ctot(e + e~ — > 7* — > hadrons) in e + e~ annihilation. In par- 
ticular at higher energies it is convenient to represent re- 
sults in terms of the cross section ratio (see Appendix A 
for more details) 



R(s) = 



ctot(e 



+ P - 



hadrons) 



a(e- 



'■y* 



fi^fi- 



(4) 



By exploiting analyticity of the irreducible hadronic vacu- 
um polarization for complex s (dispersion relation) and 
unitarity of the scattering matrix (optical theorem) it is 
possible to derive from the measured hadronic cross sec- 
tions the hadronic contribution to the photon self-energy 
II'(s). The main hadronic contributions to the shift in the 
fine structure constant is then given by |j],g| 



Aa^(s) 



3tt 



-Re 



ds 



R^(s') 



■-Ini'z. 



s f (s f — s — is) 



(5) 
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Table 1. Contribution to a£ ad = a^ d x 10 10 from exclusive 
hadronic channels and the desired accuracy for the measure- 
ment of the corresponding hadronic cross sections. 



While at large enough values of s the cross section ratio 
R(s) can be calculated in perturbative QCD, at low s one 
has to use the experimental data for R(s). A drawback 
of this strategy is the fact that theoretical uncertainties 
arc dominated by the experimental errors of the available 
e + e~ data. In (H) we have adopted the definition 



m = ^j 



Ana 2 



which is the "undressed" hadronic cross section 



T (o) 

7 had 



(s) = crhad(s) (a/a(s)) 2 



(6) 



(7) 



which agrees within 1 a with the theoretical prediction^ 



^thc 



212(190) x lfT 11 , taking (g) for , 



had 



^ ^ l - —v— ; ~ ^ , .— & VCJ , ~~ ^ ■ How- 
ever the significance of this deviation depends strongly on 
the value of a^ ad and its error Jb||l5|]. We refer to [|l6) 
for a recent review and possible implications For the near 
future a further reduction of the experimental error to a 
value of about 4 x 10~ 10 is expected, which could corro- 
borate the discovery of new physics. 

In any case the hadronic uncertainty of the theoretical 
prediction will soon be a serious obstacle for the interpre- 
tation of the expected experimental result. We therefore 
by all means need a better theoretical prediction, i.e. a bet- 
ter control of the hadronic errors. This can be by progress 
in theory as well as in more precise measurements of the 
hadronic cross sections at lower energies. 



in terms of the lowest order /i-pair production cross sec- 
tion at s> m,, . 

The procedure described is especially important for the 
precise prediction of the anomalous magnetic moment of Because of the 1/s 2 enhancement of R(s) in the inte- 
thc muon a M , to which the leading^hadronic contribution gra i (|) about 70% of the hadronic contribution to a M 

is coming from the p-io region. Not surprisingly therefore, 
the error in the prediction of a M is mainly coming from this 
low energy region. Since pion pair production at energies 
below 1 GeV is the dominant channel (see Table 0) an im- 



is given by the dispersion integral f||[|[§,@,|| 



,had 



V 3tt 



ds 



R{s) K(s) 



(8) 



lint 



proved measurement of the process e + e 



P,U3 



This integral is similar to (g|), however with a different ker- 
nel K{s), a bounded function which increases monotoni- 
cally from 0.63 at threshold (s = 4m 2 ) to 1 at s — ► oo. The 
theoretical error of a M is largely due to the uncertainty of 
the hadronic contribution || (see also [ p~0| , |ll| , |l^] ) : 



with per cent accuracy could already improve the theore- 
tical prediction of a M substantially^ 

The 7r + 7r - data are usually represented in terms of the 
pion form factor F 7r (s). The latter is related to the total 



_had 



= (697.4 ±10.5) x 10" 



10 



(9) 



Interestingly the new experimental result from the Brook- 
haven g-2 experiment p3| which reached a substantial im- 
provement in precision leads to a new world average value 



exp 



(11659202.3 ±15.1) x 10 



-10 



(10) 



1 In (TJ a 2.6 a deviation |a° xp - af c \ = 426(165) x 10" 11 
was claimed, assuming the value for ajj ad as estimated in |7J. 
Recent progress in evaluating the hadronic virtual light -by- 
light scattering contribution mjv lead to much better agree- 
ment between theory and experiment. The new result, a change 
of sign in the leading 7r -exchange contribution, was confirmed 

in g|0|3& 

2 A recent analysis of four pion production which is impor- 
tant at energies above 1 GeV was presented in J22fl 
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cross section by 



3 s 



\FA*)f 



(11) 



where /3„ = (1 — 4m^/s) 1 / 2 is the pion velocity. For the 
cross section ratio R this reads 



R^{s) 



ft 



I^ 0) WI S 



Note that 



\F^(s)f = \F n (s)\ 2 (a/a(s)) 2 



(12) 



(13) 



is the equivalent of (R) for the pion form factor. The aim of 
the present work is to discuss in some detail how to extract 
precisely the pion form factor from the experimental data. 

Present measurements are performed at the e + e~ colli- 
ders DA^NE at Frascati |j and VEPP-2M at Novosi- 
birsk @,||. While at VEPP-2M in a scan data for dif- 
ferent center of mass energies are taken, at the DA<£NE 
experiment which is running on the (f> resonance for the 
next years the radiative return due to IS photons is used 
to measure hadronic cross sections below 1.02 GeV. At 
present the experimental analysis is based on events with a 
tagged photon (2^,^,^8,29,3^]. The radiative return phe- 
nomenon also allows to measure low energy cross sections 
at the B-factories BABAR/SLAC and BELLE/KEK @. 
At higher energies R(s) measurements are performed by 
the BES Collaboration at BEPC §3. Future plans at- 



tempt to remeasure R(s) in the range M<p < E c 
(PEP-N project at SLAC). 



< M 



J/4 



In this paper in contrast to the photon tagging approach 
we focus on an inclusive treatment of all photons, includ- 
ing virtual photons which materialize into anything non- 
hadronic. This provides a cross-check of the tagged photon 
method. Furthermore, we are able to gain control over the 
theoretical error of the calculations as the fullQ 0(a 2 ) IS 
corrections are available only for the inclusive treatment 
p3{ but not for the case of a tagged photon. These correc- 
tions appear to be important since we observe large effects: 
For the pion pair invariant mass distribution da/ds' which 
is the observable measured at DA<£NE we find an effect 
of up to 15% from 0(a 2 ) IS photonic corrections and of 
up to 8% from IS pair production. The Yennie-Frautschi- 
Suura resummation of higher order soft photons and the 
leading collinear log 0(a 3 ) corrections gives us each an 
additional contribution of about half a per cent and a 
good estimate for the accuracy we can expect from our 
treatment of radiative corrections. For the tagged pho- 
ton method such an estimate seems to be more difficult 



3 The terminology used in this paper is the following: "Born 
approximation" is related to the process e + e~ — > 7* — > tt + iv~ 
without any additional photon attached to it; u O(a n ) photonic 
corrections" are obtained from the Born process by attaching 
n additional real or virtual photons to it. For the case of IS 
pair production the leading order QED corrections are already 
ofO(« 2 ). 



since, as already mentioned, the complete 0(a 2 ) correc- 
tions are still missing. Although tagging a photon has an 
advantage concerning the reduction of background, which 
is mainly coming from the processes)] e + e~ — > tt + tt~tt° 
and e + e _ — > ^ + /^~7, the theoretical uncertainty is going 
to dominate as soon as the experimental error is reduced 
to at least per cent level. The disadvantage concerning 
background reduction is partly compensated by a larger 
cross section for the inclusive method in respect to the 
tagged photon method. To obtain the pion pair invariant 
mass distribution with high accuracy the energy and mo- 
menta of the pions are measured in the drift chambers 
of the KLOE detector. The tagging of the photon is not 
necessary for this. Additionally to the IS corrections also 
FS and IFS interference corrections are considered. We 
observe large effects from FS contributions of up to more 
than 15 per cent (in the very soft and very hard photon re- 
gion) to the pion pair invariant mass distribution da/ds' . 
With the tagged photon method the FS contribution can 
be reduced by making strong cuts. However, we find that 
even for very strong cuts which reduce the cross section 
considerably the FS contribution still contributes up to a 
few per cent to the cross section. Therefore even for this 
scenario FS corrections cannot be neglected. 

One of the basic problems in calculating QED correc- 
tions to a process involving hadrons concerns the extended 
structure of the final state particles. Fortunately, the pro- 
cess we are interested in, e + e~ — > tt + tt~ , is a neutral ex- 
change channel which allows a separate consideration of 
IS radiation and FS radiation, and the latter only is cau- 
sing troubles. At long wavelength it is certainly correct to 
couple the charged pions minimally to the photon, i.e., to 
calculate the photon radiation from the pions as in scalar 
QED. In contrast hard photons couple to the quarks. Thus 
one knows the precise value of the QED contribution only 
in the two limiting cases while we are lacking a precise 
quantitative understanding of the transition region. In ad- 
dition at the p-resonance one is actually not producing a 
charged pion pair but the neutral vector-boson p , which 
further obscures a precise understanding of the radiative 
corrections. In the present paper we will first consider the 
QED corrections for point-like pions, which can be gene- 
ralized to a description of the pions by the pion form factor 
F^(s); graphically: 





plus one, two or more virtual and/or real photons attached 
in all ways to the charged lines. 



4 Note that together with the channel e + e~ — ► 7r + 7r _ 7r also 
e + e~ — + 7r + 7r~7r°* — » 7r + 7r~77 has to be subtracted as a back- 
ground. 
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Why can this procedure be trusted? There are two main 
points which convince us that the model ambiguity of the 
FS radiation cannot be too large, although we cannot give 
a solid estimate of the uncertainty. In our conclusions be- 
low we will be more concrete on this issue. The first point 
is that the FS QED corrections are ultraviolet (UV) finite 
in our case. This is in contrast, for example, to the weak 
leptonic decays of pseudo scalar mesons, where the QED 
corrections to the effective Fermi interaction depends on 
an UV cut-off, which in the SM corresponds to a large 
logarithm which probes the short distance (SD) structure 
of the hadron. There is no corresponding SD sensitivity 
in our case. This is confirmed by a recent analysis of the 
radiative correction to the pion form factor at low energies 
within the frame work of chiral perturbation theory J38| • 
In fact the correction does not depend on any chiral low 
energy parameter, which would encode an eventual SD 
ambiguity. The second important point is that the FS cor- 
rection turns out to be large (of order 10 %) in regions 
which are dominated by soft photon emission where the 
treatment of the pions as point particles is actually justi- 
fied. 

Another problem concerns the treatment of the vacuum 
polarization effect. In the theoretical prediction which is to 
be compared with the data, the photon propagator has to 
be dressed by the vacuum polarization (VP) contributions 
(for details see Appendix B): 



wvw 



n(0), 



rA/O^^ 



To extract |i 7 i u, '(s)| 2 from the experimental data one has 
to tune it by iteration in the theoretical prediction such 
that the experimentally observed event sample is repro- 
duced. Of course the appropriate cuts and detector effi- 
ciencies have to be taken into account. If one includes the 
VP effects in the theoretical prediction we obtain \F^ (s)\ 2 

or Oj^Jj while omitting them would yield l-F^s)! 2 or Chad- 
In principle, one may calculate one from the other by a 
relation like (R). The cross section ratio R(s) is only used 
as an "undressed" quantity. 

Aiming at increasing precision one has to define precisely 
which quantity we want to extract from the data. For the 
calculation of the hadronic contributions rt|) or (@) one 
must require the full one particle irreducible (lpi) photon 
self-energy "blob" which includes not only strong interac- 
tions but also the electromagnetic and the weak ones. For- 
mally the relevant quantity is the time-ordered product 
of two electromagnetic quark currents (0) which in low- 
est order perturbation theory in the SM is just a quark 
loop. While the weak interactions of quarks at low en- 
ergies are negligible the electromagnetic ones have to be 
taken into account. The leading virtual plus real inclu- 
sive photon contribution to pion pair production is about 
0.7% for s 3> 4m 2 . and increases due to the Coulomb 
interaction (resummation of the Coulomb singularity re- 
quired) when approaching the production threshold. Up 



to IFS interference which vanishes in the total cross sec- 
tion, the virtual plus real FS radiation just accounts for 
the electromagnetic interactions of the final state hadrons, 
which is "internally dressing" the "bare" hadronic pion 
form factor. Thus at the end we have to include somehow 
the FS QED corrections into the hadronic cross section. 
This means that in the photon self-energy one also has 
to include photonic corrections to the hadronic lpi blob, 
graphically: 



+ 



^/\x1iad)v/\/ 



Including the FS photon radiation into a dressed pion form 
factor F-X looks like if we don't have to bother about the 
radiation mechanism in the final state. However, it is not 
possible to distinguish between IS and FS photons on an 
event basis. Radiative corrections can only be applied in 
a clean way if we take into account the full correction at 
a given order in perturbation theory. In addition e + e~ 
pairs have to be included since they cannot be separated 
from photonic events if they are produced at small an- 
gles in respect to the beam axis. Below we will present 
and discuss the theoretical prediction for pion-pair pro- 
duction with virtual corrections and real photon emission 
in terms of a bare pion form factor F^ . We therefore ad- 
vocate the following procedure: Try to measure the pion 
pair invariant mass spectrum in a fully inclusive manner, 
counting all events 7r + 7r~, 7t + 7t~7, 7r + 7r~77, 7r + 7r~e + e~, 
7r + 7r~e + e - 7 ... as much as possible and determine the bare 
pion form factor F„ by iteration from a comparison with 
the observed spectrum to the radiatively corrected theo- 
retical prediction in terms of the bare pion form factor. 
In the theoretical prediction the full vacuum polarization 
correction has to be applied in order to undress from the 
reducible (non-lpi) effects. At the end we have to add 
the theoretical prediction for FS radiation (including full 
photon phase space). The corresponding quantities we will 

denote by F^ (s) or o^ld- To be precise 



|FW( S )| 2 HF(°>( S )| 2 (i+v(s) 



a 



(14) 



to order 0(a), where r](s) is a correction factor which will 
be discussed in Sec. 2. The corresponding O(a) contri- 
bution to the anomalous magnetic moment of the muon 
(|) is 5~<al ad = (38.6 ± 1.0) x 10 -11 , which compares to 
(46.0 ±0.5 ±9.0) x lO" 11 estimated in jn) (sec also |p|). 

One could expect that undressing from the FS radiation 
and adding it up again at the end would actually help to 
reduce the dependence of the FS radiation dressed form 
factor F-f (s) on the details of the hadronic photon ra- 
diation. We will show that this is not the case, however. 
In the radiative return scenario we are interested here, 
FS corrections depend substantially on the invariant mass 
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square s' of the pion pair and reach more than 10 % when 
s' < s (soft photons) while the FS radiation integrated 
over the photon spectrum which has to be added in order 

to obtain the 0(a) corrected pion form factor Fyf (s) is 
below 1.0%. 

The low energy determination of R(s) is complicated by 
the fact that an inclusive measurement in the usual sense 
which we know from high energy experiments is not pos- 
sible. At low energy also hadronic events have low multi- 
plicity and events can only be separated by sophisticated 
particle identification. In our case the separation of [i[i 
pairs from irn pairs is a problem which requires the ap- 
plication of cuts. However, since the /x-pair production 
cross section is theoretically very well known one may pro- 
ceed in a different way: one determines the cross section 
e + e~ — > 7r + 7r~, /i + /i~ plus any number of photons and 
e + e~ pairs and subtracts the theoretical prediction for 
e + e~ — > /i + /i~ plus any number of photons, including vir- 
tual ones materializing into e + e~ pairs, and then proceeds 
as described before. At least this could provide important 
cross checks of other ways to handle the data. 

Often experiments do not include (or only partially in- 
clude) the vacuum polarization corrections in comparing 
theory with experiment. An example is the CMD-2 mea- 
surement of the pion form factor |2J,g5|, where no VP 
corrections have been applied in determining |i^(s)| 2 Q 
The so determined form factor includes reducible contri- 
butions on the photon leg: 



This "externally dressed" form factor is not what we can 
use in the dispersion integrals. The lpi photon self-energy 
we are looking for, which at the end will be resummed to 
yield the running charge, by itself is not an observable but 
a construct which requires theoretical input besides the 
measured hadronic cross section. In fact the irreducible 
photon self-energy is obtained by undressing the vacuum 
polarization effects according to (0). A more detailed con- 
sideration of the relationship between the irreducible pho- 
ton self-energy and the experimentally measured hadron 
events will be briefly discussed in Appendix B. 

Our results are presented and discussed in the next sec- 
tion. The importance of IS and FS corrections to da/ds' 
can be seen in Fig. ||. In Fig. || da/ds' with 0(a 2 ) IS 
and 0(a) FS contributions are shown for realistic angular 
cuts. The other figures and tables are related to the inves- 
tigation of higher order photonic corrections, IS pair pro- 



duction contributions, pion mass effects, IFS interference 
corrections (da j d cos 6*) and the precision of the numerical 
calculations. A case of a tagged photon with strong kine- 
matical cuts is also briefly discussed. In Sec. 3 we consider 
the determination of | F n \ 2 by an inclusive measurement of 
the pion-pair spectrum in a radiative return scenario, like 
possible at DA^NE. Conclusions and an outlook follow 
in Sec. 4. Considerations on the experimental determina- 
tion of R(s) are devoted to Appendix A. Details about 
"undressing" physical cross sections from vacuum polari- 
zation effects are given in Appendix B. In Appendix C 
we comment on the form factor parameterization of the 
7r + 7r~ final state. 



2 Analytic and Numerical Results and Their 
Discussion 



In the Born approximation the cross section for the pro- 
cess e~(pi) + e + (p 2 ) — ► Tr~(ki) + ir + (k 2 ) is of the form 



dap 
dQ 



a^i(s) 



sm 2 e\F„(s)\< 



(15) 



5 In the final presentation of the CMD-2 data |E5J VP correc- 
tions have been applied together with the FS correction ( [14J ) 
to the "bare" cross section referred to as crJJ^f _,■> 



where is the angle between the it momentum and the 
e~ momentum, s — (jp\ + P2) 2 and (3 w (s) = yl — 4m 2 /s, 
with m^ being the pion mass. The form factor F 7T (s) en- 
codes the substructure of the pions (see Appendix C). It 
takes into account the general Tt + n~j vertex structure 
and in particular satisfies the charge normalization con- 
straint F^(0) = 1 (classical limit). In this section we will 
only consider the radiative corrections which means that 
we are considering what we denoted by F^. The vacuum 
polarization effects may be accounted for at the end via 
(|l3|). By s' = (fci + k-i) 2 we will denote the invariant mass 
square of the pion pair. 

Let us now consider the radiative corrections to the Born 
process which are related to additional virtual and real 
photons. These kind of corrections have been extensively 
studied in the literature at the one- loop level |39J, |4(3| , [4l| , |4^] 
and have been applied in the past by experiments in e + e~ 
cross section measurements. More recently radiative cor- 
rections to pion-pair production have been reconsidered 
in |4J| and were applied by the CMD-2 Collaboration for 
the determination of the pion form factor |24|,£5|. For our 
purpose, we found it necessary to redo these calculations 
for the 7r + 7r - production channel. To estimate the im- 
portance of the different QED correction contributions we 
begin with an analysis of the cross sections without kine- 
matical cuts for the total cross section a and the pion 
invariant mass distribution da/ds'. Here only IS and FS 
corrections have to be taken into account since as a con- 
sequence of charge conjugation invariance of the electro- 
magnetic interaction the IFS interference corrections do 
not contribute to these observables. 
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Initial State Corrections 






Final State Corrections 






Interference Corrections 





Fig. 1. Virtual and real 0{a) QED corrections to the pro- 
cess e + e~ — > 7T + 7r - , excluding vacuum polarization diagrams. 
The dots stand for the remaining IFS interference correction 
diagrams. 



The IS corrections include the 0(a 2 ) |33) and the leading 
log 0(a 3 ) [Q photonic corrections as well as the con- 
tri bution s from initial state fermion pair production pa , 
(see Fig. @). Among the latter only e + e~ pair 



production is numerically relevant. The FS corrections are 




,"K 



a) 



b) 



Fig. 2. Initial state fermion pair production. Diagram o) 
shows an example of a non-singlet contribution, f + f~ being a 
fermion pair which is radiated off the initial state electron or 
positron. For / = e also singlet contributions like diagram b) 
have to be taken into account. 

given to 0(a) where the pion masses are kept everywhere. 
Yennie-Frautschi-Suura resummation 44,49] was applied 
to the IS and FS soft photon contributions. We then ob- 
tain {z = s'/s): 
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2/3^ 

1+Ar 



~Sfin(s,s') 



2Li 2 
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2a M _ ^/W> 
- ( Zj /33( S ) 



2Li 2 



1-Ar 
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(37) 
(38) 



Li 2 (x) = / — log(l-y), 



Sl,2(*) 



i/^log 2 (l-y),L e = log 
27 y 



The Ofeii corrections (|23) , ( Pq ) and (|33|) are taken from ^ 
and [|§|j respectively. The dots in (§|), @, @ and 
(|32j) correspond to 0(cv 2 ) contributions which do not con- 
tain any \og(s/m 2 ) terms and can be neglected safely. 
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Fig. 3. Pion pair invariant mass distributions {da/da') with 
radiative corrections, normalized to da/ds' with only O(a) IS 
corrections. The thick line shows the case when up to 0(a ) 
IS and 0(a) FS contributions (excluding IS pair production) 
are taken into account and appropriately resummed [see ( Jlq - 
|lq)1- The thin solid line shows the same but this time with- 
out resummation. The dotted line corresponds to the 0(a) FS 
corrections (together with O(a) IS corrections). For the long- 
dashed and the dot-dashed lines only the resummed IS 0(a 2 ) 
and the resummed IS O(a) radiative corrections are taken into 
account, respectively. 



Fig. H shows the pion pair invariant mass distributions 
da/ds' with radiative corrections normalized to da/ds' 
with only O(a) IS corrections (y/s = 1.02 GeV). In Table | 
the contribution from IS 0(a 2 ) and FS O(a) photonic cor- 
rections are shown for different center of mass energies. 

The following points can be recognized: 




Fig. 4. Pion pair invariant mass distributions (da/ds') for 
different center of mass energies y/s = 0.6,0.76,0.9,1.0,1.02 
GeV. The solid lines stands for the "complete" cross section, 
including 0(a 2 ) IS and 0{a) FS corrections [see (0-0)]. The 
dotted lines give the results when the 0(a 2 ) IS corrections are 
neglected. The dot-dashed lines correspond to the case when 
the 0(a) FS contribution is neglected. 



The FS corrections (dotted line) are quite large, espe- 
cially in the region of soft photons as well as for very 
hard photons; 

The 0(a 2 ) IS effects are considerable; 
The FS and IS contributions compensate each other 
significantly for large s' . 



Fig. H shows da/ds' for different center of mass energies 
*Js. Going to smaller center of mass energies the 0(a 2 ) 
IS corrections become smaller and smaller. On the other 
hand the FS contributions remain considerably large. In- 
terestingly, for the 4> resonance energy (y/s — 1.02 GeV) 
both the 0(a 2 ) IS and the 0(a) FS contributions are 
large. Quantitatively this is shown in Table 0. The re- 
summation of the 0(a 2 ) IS soft photon logarithms [see 
(|l7|)1 gives a contribution smaller than 5 per mill for s' 
below the p resonance peak and smaller than 3 per mill 
above it. The resummation of FS soft photon logarithms 
[see (fig)] changes the complete results only slightly (less 
than 0.5 per mill). To reduce the theoretical error to a 
few per mill one also has to include the contributions from 
initial state e + e~ pair production |33J, [35| , p6|,|37| , given in 
@ and <M). In Table § and g the Oia 2 ) and leading 



) pair production contributions to da/ds' for differ- 
ent hadronic energies are presented. What is remarkable 
is the very large singlet contribution [see Fig. |]b)] in the 
region of low hadronic energies which amounts about 8 
per cent for yfs' = 0.3 GeV. Since these effects are related 
to e + e~ pairs which are mainly emitted collinearly to the 
beam axis they escape detection and therefore have to be 
included into the data analysis. Hence when unfolding the 
data from radiative corrections also these effects have to 
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Vs 1 [GeV] 


0(a 2 ) IS 


0(a) FS 




contribution 


contribution 


0.3 


4.3 


11.5 


0.4 


4.4 


4.3 


0.5 


4.0 


3.2 


0.6 


3.4 


2.0 


0.7 


2.2 


0.9 


0.76 


1.2 


0.7 


0.8 


0.2 


1.2 


0.9 


- 3.6 


5.5 


0.95 


- 6.9 


10.1 


1.0 


- 15.3 


16.6 



Table 2. Contribution of 0(a 2 ) IS and O(a) FS corrections 
to da /ds' (in %), Vi = 1.02 GeV. 



Vs 1 [GeV] 


0(q 3 ) IS 




contribution 


0.3 


3.9 


0.4 


4.3 


0.5 


4.2 


0.6 


3.8 


0.7 


3.0 


0.76 


2.4 


0.8 


1.8 


0.9 


0.3 


1.0 


0.6 



Table 5. O(o 3 ) 
(in per mill). 



leading log IS photon contribution to da/ds' 



v^ 7 [GeV] 


0(a 2 ) IS pp 


Singlet 




contribution 


contribution 


0.3 


79.1 


74.9 


0.4 


36.3 


31.9 


0.5 


16.6 


12.2 


0.6 


8.3 


4.0 


0.7 


4.8 


0.76 


0.76 


3.9 


- 0.01 


0.8 


3.4 


- 0.24 


0.9 


2.7 


- 0.27 


1.0 


1.2 


0.06 



Table 3. 0(a 2 ) contribution from IS pair production to da/ds' 
(in per mill). In the second column only the singlet contribution 
(including singlet-non-singlet interference) is shown. 



Vs 1 [GeV] 


0(a 3 ) IS pp 


Singlet 




contribution 


contribution 


0.3 


- 0.87 


- 1.27 


0.4 


- 0.45 


- 0.82 


0.5 


- 0.18 


- 0.48 


0.6 


- 0.07 


- 0.27 


0.7 


- 0.09 


- 0.15 


0.76 


- 0.15 


- 0.10 


0.8 


- 0.21 


- 0.07 


0.9 


- 0.43 


- 0.03 


1.0 


- 0.72 


- 0.001 



Table 4. 0(a 3 ) IS pair production contribution to da/ds' (in 
per mill). 



be subtracted. The leading contribution from 0(a 3 ) pair 
production appears to be less than 1 per mill which gives 
us a good estimate about the precision we can expect. 

We also take into accout the leading log O(ar) IS photon 
correction J34|, which is given in (|4|) and (|2q). The con- 
tribution can be of the order of 4 per mill for hadronic 
energies below the p resonance peak, as shown in Table pi 

The total cross section a(s) can be obtained by carrying 
out the s' integration in (|16|) numerically. For a(s) the 



0(a 2 ) IS corrections are not as important as for da/ds' 
(they account for at most 1 %, at y/s — 1.02 GeV 2 per 
mill). Neglecting the 0(a 2 ) corrections, cr(s) can then be 
written in the following simple form: 



cr(s) = cr [1 + S ini (A) + 6 fm (A)} 

3-2y/sA 

ds' <T (s')p ini (s,s') 



4 m 2 

o-q(s) 



(39) 



-2V5A 



4/H^ 



ds' Pf in (s,s') , 



with 



S ini (A) = log(— =j B e (s) 
'2A 



EV-bS(l), 
ini 



SLT jyL> (s) 



(40) 



Sfin(A) = log 
ini \Sj S ) 



V+S(l), 



B„(s,s')+6 v f ^>(s), (41) 



fff(i) 



C W) + 



B e (s) 



Pfvn{s,s') = - 



1 1 — z 



(42) 
(43) 



where A is the soft photon cut off energy which drops out 
in the sum (B9h. 



The total cross section is plotted in Fig. 0. In Table @ 
the 0(a) FS corrections to a(s) for different center of 
mass energies are shown. Although it can be hardly rec- 
ognized directly from the figure, the FS contributions are 
not marginal for energies below the p resonance peak. 

Taking the high energy limit in ( p9[ ) provides a good cross 
check for the FS correction results. Carrying out the s' 
integration for s — > oo and adding this result to the high 
energy virtual and soft photon FS corrections leads to an 
expression in which the s dependence drops out. This has 
to be the case according to the Kinoshita-Lee-Nauenberg 
theorem J46|,[47j which requires the collinear logarithms to 
cancel. In addition all terms proportional to ir 2 drop out. 
Defining 



V(s) 



S fin {A) + 



-2^/sA 



iin'z 



ds' pf in (s,s') 



, (44) 
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- - without FS contribution 

Born only 

— complete 



XI 

lE 

C 600 - 
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0.2 




this limit we observe: 



0.6 0.S 

sqrt(s) [GeV] 



1.2 



Fig. 5. Total cross section a(s) as a function of the center 
of mass energy. The solid line corresponds to cr(s) as given 
in (|39J). The dotted line corresponds to the Born cross section. 
The dot-dashed line corresponds to the Born cross section with 
only 0(a) IS corrections. 



v^ [GeV] 


O(a) FS 




contribution 


0.3 


3.6 


0.4 


1.2 


0.5 


0.9 


0.6 


0.9 


0.76 


0.7 


0.9 


0.4 


1.02 


0.3 



Table 6. Contribution of O(a) FS corrections to the total 
cross section (in %). 



one can write the total cross section with only 0(a) FS 
corrections in the following compact way: 



o-fin(s) = 1 + 77(5)- <t (s) 

IT - 



The function r)(s) is given by @,gg,|§ 

,(.s) = ^{4Li 2 (T^)+2Li 2 (-T 



(45) 



-ft 
+ft 



-3 log 



(wk) log(£fe)-21ogG0,r)bg( 



1+ftr 
1-ftr 



-3 log (jdW 



4bg(Ar) 



+^[|(l + ^) 2 -2]log(i±f) + 



2 01 



and provides a good measure for the dependence of the 
observables on the pion mass. Neglecting the pion mass 
is obviously equivalent to taking the high energy limit. In 



i](s —» 00) = 3 



(46) 



Our result in ( 48) agrees with the result obtained by Schwin- 
ger Wq| but disagrees with that in J43| for which in this 



limit the terms 



do not drop out. In Fig. |fj r](s) 



is 

plotted as a function of the center of mass energy. It can 
be realized that for energies below 1 GeV the pion mass 
leads to a considerable enhancement of the FS corrections. 
Regarding the desired precision, ignoring the pion mass 
would therefore lead to wrong results. 

Close to threshold for pion pair production (s ~ 4m 2 ) the 
Coulomb forces between the two final state pions play an 
important role. In this limit the factor rj(s) becomes sin- 
gular [r](s) — > 7r 2 /2/3 7r ] which means that the 0(a) result 
for the FS correction cannot be trusted anymore. Since 
these singularities are known to all orders of perturbation 
theory one can resum these contribution, which leads to 
an exponentiation |4q] : 

/ a ira\ ira 



1 — exp 






(47) 



Above a center of mass energy of y/a = 0.3 GeV the ex- 
ponentiated correction to the Born cross section deviates 
from the non-exponentiated correction less than 1 %. 




1.0 1.5 

sqrt(s) [GeV] 



2.0 



Fig. 6. The FS correction factor 77(g) as a function of the 
center of mass energy yfs [see (kjH4q)]. 



We now consider the IFS interference corrections (see Fig. [j]) 
which modify the angular distribution. To 0(a) we can 
write: 



da 
dQ 



dap 
dO 



1 + 8{A) 



(48) 
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dah \ 



(A) 



where the correction factor 8 is the sum of the virtual plus 
soft photon IS, FS and IFS interference correction factors: 



5(A) = 5 ini (A) + 6 fin (A) + 8 in t(A) 



(49) 



dat/dfi is the hard photon contribution which is calcu- 
lated numerically. 8mi(A) and 5f in (A) are given in ( p0[ ) 
and (ffl[). 5i n t(A) can be written in the following, compact 
way: 



Oint = < log 

n \ 

1 s 2 



—u 



mt 



-t + ml 

u — t 



M^ 



4at — m: 
u 4- * r 
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c\(u) + Cl(t) + C e + Ctt 
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(50) 
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Li 2 
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Li 2 

' 771-772 

7/2 - b x 



772 - a 

7/2 - b a 



772-771 



with (a; = u,t, Ki — K%(x)) 



1 



ci(x) = 2I 31o S 



^--f 



!*(_!*)_„,(_£)*" J 



2Li 2 I - 

1 + /3 2 - 

+Li 2 



7T" 

T 



2Li 2 



1+/3. 



2Li 2 



1-/3. 



1 + Ar 

1-/3. 



Li 2 



1-/3. 

'l + Ar 



lo s I — ) lo s 



'»: 



1-/3. 



Ce 4 log2 (i) + T 



^(S) = [A*(«l) - /f M - /f («*) + /f(« 3 ) 

- jT («1> K 2) - /f («4, «l) + /£ («4, «2) 
+ ft («2, «l) + /l(«l, K4) + /|(«2, K4) 

- / 5 x (k 3 , «l) + /£ («3, «2) " /s («8, «a) 

- /? («1, «3) - /e ( K 2, «&) + / 6 X (K4, «*)] 



/f('7) = 5log 2 [& !C -77]-ilog 2 [a-7 7 ], 
/ 2 x ('7) = 5log 2 [77-o]-ilog 2 [77-6 :c ], 



/f (771,772) = -Li 2 
+ Li 2 f- 



(6a - a) (771 -772) 

(&x - 77i)(a - 772) 
fo T — a 



a- 772 



Li 2 



&t - 



771 



/3.(s), 6 x = /3 e (s) + 2 



«i,2 = «i,2(a;) = -1 + 



1 



V— sir 



2 2 
—x + m„ — m„ 



«3,4 = K3,4(a;) = 1 + 



± ^J\(x,ml,ml) 
1 



V-sa; 



2 2 

—x + m„ — m„ 



± \J\{x,m 2 e ,ml) 
X(x, y, z) = z 2 + y 2 + z 2 — 2xy — 2xz — 2yz. 



From ( p0| ) it can be seen immediately that 5i n t is anti- 
symmetric, thus it changes sign under the exchange £ <-> w 
[t(u) = (pi — fc2(i)) 2 ]- This is actually required by charge 
conjugation invariance. 



Table ffl shows the IFS interference contribution to the pion 
angular distribution (Fig. [7]). One can recognize, that with 
an angular cut between the pion momentum and the beam 
axis of 20° < O < 160° this is not bigger than 5.6%. 
The importance of the interference contribution can be 
enhanced by tagging the photon and imposing a strong 
cut on the ang le between the photon momentum and the 
beam axis [Eq| (see Fig. g) . This seems to be the only way 
to tackle the imaginary part of the pion form factor. 

The results presented so far have been obtained by the 
dedicated Fortran program A^pwDiTE . It generates cross 
sections with the option of kinematical cuts |5(J as needed 
by experiment. As shown before, the 0(a 2 ) IS (photonic 
and IS pair production) contributions to da/ds' are con- 
siderable and even 0(a 3 ) leading log contributions should 
be included. Analytical formulae for the full 0(a 2 ) IS cor- 
rections with cuts have not been calculated so far. At this 
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■§ 



CD 

X 


o 

T3 


"O 




0.5 

cos(9) 

Fig. 7. 7T _ angular distribution for y/s = 1.02 GeV. The solid 
line, corresponding to the complete O(a) corrections, is not 
symmetric as a consequence of the IFS interference corrections. 
Both the tree level distribution and the distribution with only 
IS and FS corrections are symmetric. 



cos© 


d<TB + IS + FS 1 drrtot 
d cos & Id cos & 


IFS interference 


-1.0 


0.99/ 1.46 


47.5 


-0.99 


2.95/3.33 


12.9 


-0.94 


11.12/11.74 


5.6 


-0.6 


44.03/44.97 


2.1 


-0.2 


59.91/60.3 


0.6 


0. 


61.76/ 61.76 


0.0 


0.2 


59.91/ 59.53 


-0.6 


0.6 


44.03/43.09 


-2.1 


0.94 


11.12/10.49 


-5.6 


0.99 


2.95/2.56 


-13.2 


1.0 


0.99/ 0.52 


-47.5 



Table 7. Contribution of the interference terms (in %) to the 
differential cross section (corresponding to the solid and dotted 
line in Fig. |7|). 



stage we therefore rely on the complete results without 
cuts ckr( comp ^ I 'ds' as given in ( |l6| ) with the following ap- 
proximationH: 



/ dfJ (compl) 

I IP 



^ ds' 



(da^ 



( d<r<.") \ 
V ds ' J, 



V ds' 



(51) 



da^/ds' is the differential cross section to O(a). See 
Fig. for an example. In the limit s' —> s the above ap- 
proximation is exact since then the radiated photons are 
soft. In principle we can expect that away from this limit 



6 For an exact treatment we would need the analytic expres- 
sion for the angular distribution at 0(a 2 ). 



-Q 
a 



c 
-a 

"D 




Fig. 8. Lowest order n" angular distribution for the case of a 
tagged photon. The angular cut between the photon momen- 
tum and the beam axis is chosen such that only photons in the 
angular range 60° < <9 7 < 120° are detected. The difference 
between the solid and dotted line is due to an additional cut 
between the tagged photon and the pions 7° < <9 77r < 173° 
(solid line). This cut was also applied to the remaining dashed 
and dot-dashed line. The curves correspond to different values 
for the minimal photon energy A. The solid and the dotted line 
correspond to A = 0.01 GeV, the dashed line to A = 0.02 GeV 
and the dot-dashed line to A — 0.03 GeV. 




0.3 



0.5 



0.7 



0.9 
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Fig. 9. Pion pair invariant mass distribution with an angular 
cut 30° < < 150° between the n momenta and the beam 
axis, Jl = 1.02 GeV. 



the situation is different since the contribution from a sec- 
ond hard photon could distort the angular distribution of 
the pions. The distortion however remains below 1 per mill 
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s' 


Set A 


Set B 


set A — set B 


all 


noFS 


all 


noFS 


all 


noFS 


0.8 


11.994 


9.981 


18.231 


16.121 


6.237 


6.140 


0.85 


12.252 


9.494 


18.201 


15.313 


5.949 


5.819 


0.9 


14.212 


10.168 


20.615 


16.384 


6.403 


6.216 



Table 8. da/ds' in [nb/GeV 2 ], for some values of s' . The FS 
contribution for a strong cut scenario (Set A-Set B) is shown. It 
is 1.6 %, 2.2 %, 2.9 % for s' = 0.8, 0.85, 0.9 GeV 2 , respectively. 

for s' > 0.3 GeV 2 and an angular cut between the pion 
momenta and the beam axis of less than 30 degreesn. 

In the case of a tagged photon the FS corrections can be 
reduced by applying strong cuts between the photon and 
the final state particles. See Fig. [n] for such a strong cut 
scenario at the </> peak. It can be seen that the strong 
cuts reduce the FS contribution considerably. However, 
as shown in Table g, the FS contribution still amounts 
up to a few per cent. Although the presented results are 
based on an 0(a) calculation (a similar approximation as 
the one given in (pW is not possible) it is highly unprob- 
able that the situation will improve if 0(a 2 ) corrections 
are included. Finally a few remarks about the A^pwDiTE 



> 

a 
■ — 

Xl 







"D 




s' [GeV 2 ] 

Fig. 10. Pion pair invariant mass distribution da/ds' for the 
case of a tagged photon. Set A corresponds to a 7° angular 
cut between the photon momentum and the beam axis and a 
30° cut between the tv momenta and the beam axis. For set 
B the pion cuts are the same but the photon cut is now 20°. 
Taking the difference (Set A-Set B) the photon is restricted to a 
region well separated from the pion momenta. The solid lines 
correspond to the complete cross section (Born plus IS and 
FS bremsstrahlung) , for the dotted lines FS bremsstrahlung is 
neglected. 



A [GeV] 


a [nb] 


Sa [nb] 


0.1 


94.907 


0.0095 


0.01 


99.123 


0.0104 


0.001 


99.394 


0.0129 


0.0001 


99.420 


0.0157 


10~ 5 


99.422 


0.0210 


10 -6 


99.422 


0.0255 


io- 7 


99.422 


0.0303 


io- 8 


99.421 


0.0357 


1Q -9 


99.422 


0.0418 


io- 10 


99.421 


0.0493 



Table 9. Cut-off dependence of the total cross section a ob- 
tained from 4-dimensional numerical integration, ^/s — 1.02 
GeV. Sa is the absolute numerical error to a. 



A [GeV] 


a [nb] 


Sa [nb] 


0.1 


94.909344406421 


2-10 _y 


0.01 


99.126309344279 


2- 10 -11 


0.001 


99.396403660854 


2 • IO" 9 


0.0001 


99.422466900996 


3 


IO" 9 


10" 5 


99.425064117054 


6 


IO" 9 


IO" 6 


99.425323747942 


7 


IO" 9 


IO" 7 


99.425349708976 


7 


IO" 9 


IO -8 


99.425352318987 


1 


10~ 8 


IO" 9 


99.425352327085 


1 


IO -8 


io- 10 


99.425352168781 


1 


IO -8 



Table 10. Cut-off dependence of the total cross section a ob- 
tained from 1-dimensional numerical integration, y/s — 1.02 
GeV. Sa is the absolute error to a 



numerically to obtain the total cross section without cuts 
(see Table ||). This is then compared to the total cross 
section obtained from ( p9| ) by one-dimensional integration 
(Tabl ejl0| ). We observe excellent agreement. Table nfl and 
Table |m addition show the total cross section as a func- 
tion of the soft photon energy cut-off A. For values of 
A < 10^ 4 GeV we get stable cut-independent results. 



3 The Pion Form Factor from Radiative 
Return 



The pion-pair invariant mass spectrum in scalar QED 
may be written in the form 



da / da 
d? = Ids 7 



da 
ds 7 



da 



(52) 



Considering only the 0(a) contribution we can write 



program. To check the numerical accuracy, the four di- 
mensional phase space integration has been carried out 

7 We thank S. Jadach for help in checking this with a dedi- 
cated MC program based on fell. 



da 



ds' 



— =N ini ( S ,s')\F n ( S ')\' x 



dcos0 7 dcosO 



Ei-^i 



point 1 2 

'"■■■'<■■ / , \ J v l ini ' 

A 
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da 

Is 1 



N int (s,s') 2Re F„(s')F;{s) x 



cuts 

da 
ds 



dcos<9 7 dcosO^- &f>«- ^2Mf£ nt Mll 



point 



N &n (s,s')\F n (s)\ 2 x 



fin 



/ dcOS(9 7 dcOS&K- d<fin- /J l^fin 
./cuts A 



point 1 2 



l^7r( s ')| 2 drops below \F 7r (M] ) )\ 2 and a precise and model 
independent determination of F w becomes more difficult. 
Note that because of the 1/s 2 enhancement in the disper- 
sion integral (|8|) the low energy tail is not unimportant as 
a contribution to a^ ad . 



4 Final remarks and Outlook 



where the Ni(s,s')'s are appropriate normalization fac- 
tors. 0„— is the ir~ production angle and <9 7 the angle 
between the emitted photon and the n~ in the center of 
mass system of the n + ir~ pair. Cuts in the laboratory 
system may be implemented easiest by first performing 
a boost from the center of mass system of the pion pair 
to the laboratory system. If the integration over O n - is 
performed with symmetric cuts in the acceptance angles 
n ± in the laboratory frame, the 0(a) interference term 
drops out due to C-invariance and we are left with the IS 
and FS terms onlyn. Photons are assumed to be treated 
fully inclusively, i.e., we integrate over the complete pho- 
ton phase space and thus obtain: 



da 
ds 1 



\F„(a 



'\|2 



sym— cut 



da\ 



point 



ds'J i 



+ \F«(a)f 



da\ 



ini, sym — cut 
point 



^ uo / fin, sym-cut 

and hence we may resolve for the pion form factor as 



\F„(s 



l\\2 



fda_\P° lnt 

\ ds' J ini, sym— cut 



\FAs)f 



da_\ 

^7 fin, 



' da 
ds 7 



point 



sym— cut 



(53) 



sym— cut 



This is a remarkable equation since it tells us that the 
inclusive pion-pair invariant mass spectrum allows us to 
get the pion form factor unfolded from photon radiation 
directly as for fixed s and a given s' the photon energy 
is determined. The point cross sections are assumed to be 
given by theory and da /ds' is the observed experimen- 
tal pion-pair spectral function. In spite of the fact that 
both terms on the r.h.s. of ( |53| ) are of 0(a) the second 
one can be treated as a correction because the IS radi- 
ation dominates in comparison to the FS radiation. We 
observe that in the determination of ^(s')! 2 via the ra- 
diative return mechanism the to be subtracted FS radia- 
tion only depends on |F^(s)| 2 at the fixed energy s = M?. 
Note that we also benefit from the fact that \F 7T (M?)\ 2 is 

small in comparison to li^s')! 2 in the most relevant re- 
gion around the p-peak. Below about 600 McV, however, 



Note that C-invariance does not forbid a C-symmetric 
0(a 2 ) interference contribution. However we can expect such 
a contribution to be negligible. 



Experimental data on pion pair production in low energy 
e + e" collisions of percent level accuracy are availbale now 
from Novosibirsk and will be available soon from Frascati. 
That is why theoretical calculations of at least an accuracy 
of the same order are needed. In this paper we presented 
and discussed analytic and numerical results which should 
allow us to reach the desired accuracy for the appropriate 
observables. We advocated to look at the 7r + 7r~ invariant 
mass spectrum in an inclusive way for what concerns the 
accompanying photon radiation. We observe that 0(a) 
massive FS corrections as well as 0(a 2 ) IS photonic and 
e + e~-pair production corrections have to be taken into 
account. Also the resummation of higher order soft pho- 
ton logarithms and leading 0(a 3 ) IS photonic and pair 
production contributions may be necessary. 

Another background which should be estimated more care- 
fully is pion pair production via the two photon process 

77 — > 7r + 7r~ 



541 (see Fig. [II]). These events have 
a different topology, typically the pion pair appears to 
be boosted in beam direction, and may be eliminated by 
appropriate event selection. At the level of total cross sec- 
tions e + e _ — ► e + e _ 7*7* — ■> e + e~7r + 7r~ is at least an order 
of magnitude smaller than the leading pion pair produc- 
tion mechanism at energies below the <p mass. 
Supposing that the two-photon 7r + 7r~ production is, or 



53 



7"V-7T- 



+ 

- e~ e~ 



Fig. 11. The two-photon pion pair production mechanism. 

can be made, sufficiently suppressed, and under the con- 
dition that pion-pair acceptance cuts are applied in a C- 
symmctric way and hence the IFS correction drops out, 
the inclusive pion-pair distribution da / ds' is of the form 
© 

— = a (s') pini(s, s') + CT (s) pfi n (s, s') , (54) 

which we may solve for ao(s') (alternative form of (p 



ao(s') = 



da 



Pini(s,s') [ds 



7 - C (s) Pfin(s,s') 



(55) 



At DA^NE s is fixed at s = Ml and hence the FS ra- 
diation factor multiplies the fixed pion-pair cross-section 
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<tq(s = M?) at the <p. The FS subtraction term in (55) is 
an at most 10 % correction of the first and leading term for 
0.3 GcV < \~P < 0.95 GeV (in the p resonance region the 
contribution is of the order of 1 %), although both terms 
are formally of the same order 0(a). 

Such a measurement should be complementary to the pho- 
ton tagging methodn, which is not yet as well under con- 
trol as the inclusive pion mass spectrum. Since the pro- 
cess e + e~ — > /i + /i~ is theoretically very well under con- 
trol but the separation of ir + n~ and fi + p~ states is quite 
non-trivial, experimentally one actually should perform 
an inclusive measurement also with respect to muon pair 
production and then subtract the theoretical /i + /i~ cross 
section. At least this could provide an important cross 
check for the particle identification procedure. 

Apart from the fact that it would be desirable to have 
available a full 0(a 2 ) calculation for the differential cross 
section, the main limitation of our approach lies in apply- 
ing scalar QED to the pions generalized to an arbitrary 
pion form factor up to non-factorizing 0(m^/s) effects. 

We would like to stress once more that there are strong in- 
dications that the treatment of point-like pions together 
with its generalization to extended pions modeled by a 
form factor provides a reliable framework for extracting 
the pion form factor from the data. The sensitivity to the 
quark structure is minimized for the relevant observables 
by the fact that the QED radiative corrections are ul- 
traviolet finite and hence no large renormalization group 
log's show up. Furthermore, the region s' < s exhibiting 
large FS corrections corresponds to the soft photon regime 
where our generalized scalar QED treatment of the pho- 
tonic corrections is reliable. However, the fact that the 
corrections which could be sensitive to the hadronic com- 
positeness are small in the region where hard photons are 
involved does not mean that uncertainties are small at 
low s' . The reason is that for small s' the emitted pho- 
tons are hard and therefore can probe the substructure of 
the pions. One therefore can question the applicability of 
scalar QED when treating FS radiation in this region. At 
the same time it is the region where l-F^s')! 2 drops below 
|i^(M?)| 2 which enhances the FS contribution in (|53|). 
The uncertainty in the FS correction term carries over to 
the extracted form factor. 

Let us mention that the fact that we have to include FS 
corrections according to ( |l4| ) does not reduce the sensiti- 
vity to the details of the emission of photons by hadrons, 
because the FS correction one has to subtract (see (p3|)) 
is different from what one has to add at the end. The 
first reflects the photon spectrum locally, the second is an 
integral over the photon phase space. 

As a crude estimate of the uncertainty related to the pion 
substructure we replace the pions by fermions of the same 



charge and mass F]. Hence in (pq) in stead of pfi n we take 
the fermion final state radiator function 



pLM 



«l + 7 /W) 



s' &(«) s + 2m 2 { 0„(s') 



x log 

- 1- 



1 + flrQO 
1 - Ar(s') 

4s'm 2 



1 -Ami- 



2m 



2 1 



(56) 



Results are shown in Figs. (0,0). In the soft photon 
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Fig. 12. Worst case estimate of uncertainty in FS radiation 
due to non point-like structure of the pion: Fermionic vs. scalar 
radiator function. 
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Fig. 13. Influence of FS radiation on the extracted pion form 
factor (worst case estimate). 

region we have Pt in {s' < s) ~ pfin(s' ~ s) which reflects 
the correct long range behavior. For the extraction of the 
pion form factor we observe deviations of the fermionic 
from the scalar approach of less than 0.1 % for energies 
above 560 MeV. For y/s* > 420 MeV the deviation is less 



For recent progress see |55l 



10 We cannot just replace the pions by the quarks produced 
in first place because the wrong net charge would not allow to 
match the proper long distance limit. 
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than 1 %. At lower energies the difference between both 
approaches becomes larger since the radiated photons be- 
come harder: at yfs' = 360 MeV we observe a deviation 
of 2%, at Vs 7 = 300 MeV of 6.5% which is of the same 
order as the complete FS contribution in this region. Con- 
cerning the determination of a^ ad we obtain a difference 
between the fermionic and the scalar approach of about 
2(7) per mill if we restrict the analysis to a region where 
V7 > 420 (300) MeV (see Fig. (@). 




Fig. 14. FS radiation uncertainty of ajj ad as a function as a 
function of the lowest energy data point accepted (worst case 
estimate) . 



The guesstimate looks reasonable because the such ob- 
tained uncertainty goes to zero in the classical limit (s' — > 
s) and becomes of the order of the FS radiation itself in 
the hard photon limit. Note that the increasing uncer- 
tainty for low energies vV here is a consequence of the 
radiative return method since in this region the emitted 
photons are necessarily hard. 

The error due to the missing FS 0(a 2 ) and IS 0(a 3 ) cor- 
rections (including initial state pair production contribu- 
tions) is estimated to be not more than f per mill, re- 
spectively. Concerning the QED corrections we therefore 
estimate the accuracy to be at the 2 per mill level. On top 
of the perturbative uncertainty we have to take into ac- 
count the hadronic uncertainty discussed in the previous 
paragraph. 
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Appendix A: Experimental determination 
of R(s) 

As mentioned in the introduction the hadronic cross sec- 
tions are conveniently represented in terms of the cross 
section ratio 



R(s) 



__ Ohad(s) 



r (0) 
'had 



(*) 



(0) 



(*) 



(57) 



are not di- 



The "physical" cross sections Chad and c w 
rectly observable but are the result of the usual unfolding 
from real and virtual photon radiation. The "undressed" 
cross sections a\ (s) are related to the physical ones by 

<j\ (s) = ai(s) (a/a(s)) 2 [||. Obviously the effective cou- 
pling a(s) entering the physical cross sections drops out 
from the cross section ratio and hence may be replaced by 
its low energy value a. Here we briefly discuss how experi- 
ments determine R(s) and what the problems are thereby. 
By the relation ( [f2| ) our comments apply to the pion form 
factor as well. 

A direct measurement of the ratio of the physical cross 
sections ffhad and c MM has the advantage that certain un- 
wanted effects drop out from the ratio. This in particular 
concerns the normalization and its uncertainties but also 
the vacuum polarization effects. What still has to be cor- 
rected for is phase space of ct^ in the threshold region 
and the difference in final state radiation. 

Experiments up to now do not actually determine the ratio 
of the physical cross sections Chad and er MM . By the usual 
limitations in statistics and the fact that o^ drops like 
1/s not far above threshold it would not be an optimal 
strategy to do so. 

In practice at first the integrated luminosity for each mea- 
surement must be determined from the measurement of a 
reference process like Bhabha scattering, typically. Thus 
experiments in fact determine 



D , v Nhad (f + 8rc) ct : 
R{S) = N 1 ~ 

1 v norm c 



norm 



» 



a 



/<■/<■ 



o(s) 



(58) 



from the ratio of the number of observed hadronic events 
Miad to the number of observed normalizing events N noim . 
The correction <5rc incorporates all radiative corrections 
to the hadron production process, e is the efficiency - ac- 
ceptance product of the hadronic events and cr n orm(s) is 
the physical cross section for the normalizing events in- 
cluding all radiative corrections integrated over the accep- 
tance used for the luminosity measurement and tr MM . q(s) = 
47ra 2 /3s. 

This shows that the determination of R(s) depends a lot 
on the theoretical state of the art calculations used to ana- 
lyze the data. Unaccounted radiative corrections, or sim- 
plifications often made in view of other sources of uncer- 
tainties, usually contribute substantially to the systematic 
error of a measurement. 
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As mentioned above life would simplify a lot if one would 
have high enough statistics which would allow us to apply 
the definition (|57]) directly to the data. This method would 
be suitable for a precise determination of the low energy 
tail of the pion form factor where <7 MM has not yet dropped 
too much. In this region one actually should redefine R(s) 
by (see (§)) 



R(s) 



f \ —DO 

_ CThad(s) G^ 



Born 



a 



.(*) 



47TQ 2 



O-had(s) 



,(*) 



1 - Am 2 /s (1 + 2m 1 / s) (59) 



in order not to introduce fake phase space effects which 
have nothing to do with the hadron cross section which 
R is supposed to represent. The dispersion integral repre- 
sentations of Z\ev had and aj^ ad in terms of R(s) otherwise 
would have to be modified appropriately. 

Except from the threshold region one has to rely on the 
much more involved procedure described above. 



propagator which yields an on-shell version of the effec- 
tive electromagnetic coupling. In the latter approach also 
non-logarithmic contributions are resummed and in gene- 
ral this makes sense only if these contributions are small 
enough such that it does not matter whether one takes 
them into account in resummed or in perturbatively ex- 
panded form. This usually is the case at high energies 
where the log's are large and clearly dominate. The dif- 
ference between the RG and the Dyson resummation ap- 
proach is less problematic in the Euclidean (space-like or 
t-channel) and it is common practice to work with the 
space-like effective charge and to take into account the 
terms specific to the time-like region separately. In pertur- 
bation theory the difference is given by the i7r-terms from 
logarithms with negative arguments: — log(— q 2 /m 2 ) = 
^(\og(q 2 /m 2 ) — in). Since n'^s) is complex at s > 4m 2 ± 
(or above m 2 when 7r°7 production is included) one could 
consider a complex a(s) (see (g) and @)) via the shiftp] 



Aa(s) = Aa lcp (s) + Aa^(s 



(62) 



with 



Appendix B: Vacuum polarization 



Here we comment on problems related to the treatment of 
the vacuum polarization corrections at low time-like mo- 
menta. They must be included in order to avoid unneces- 
sary additional systematic errors which could obscure the 
interpretation of the experimental results. In principle, for 
processes like pion pair or muon pair production this can 
be easily accomplished, namely by applying (R) to the 
physical cross section obtained by unfolding from the other 
QED corrections (as presented in this paper). However, 
in the reference process needed for the normalization the 
situation in general is much more involved. For example, 
if wide angle Bhabha scattering is applied for the lumino- 
sity monitoring, there are different scales involved due to 
the mixed s- and t-channel dependences. Thus the depen- 
dence of measurements like R(s), or equivalently l-F^s)! 2 , 
on vacuum polarization effects is rather complicated as the 
effective fine structure constant enters in various places 
with different scales. Vacuum polarization effects if not 
accounted for properly in the data analysis are thus hard 
to reconcile at a later stage. 

There is another problem: in applying (uh formally a(s) is 
required at low energies in the time-like region. However, 
particularly in the resonance regions, this is a strongly 
varying function defined by the principal value (PV) inte- 
gral d|) and (nj). In regions where R(s) is given by data, 
the PV integral is quite ill-defined and one would have 
to model or smooth the data before integration. In ad- 
dition, the running coupling a(s) has to be seen in the 
spirit of the renormalization group (RG), which in first 
place is a systematic summation of the leading-log's, the 
next-to-leading-log's, and so on. The definition via (||) and 
(|j]), adopted commonly for the effective fine structure con- 
stant, corresponds to the Dyson summation of the photon 



Aai ep (s)= Y, (^a i i ) {s) + ^a l t\s) 



t=e,fi,T 



/>oo 

I \ as / / 

Aa had (s) = -— Re / ds —— 

3tt J Am 2 s'(s' 



R(s') 



ie) 



i^R(s) 



(63) 



The imaginary part of Z\a had is directly proportional to 
R(s). Thus in perturbative QCD R(s) ~ N C J2Q} (1 + 
0(0,5/71-)) with N c the color factor. Non-perturbative con- 
tributions to R from resonances may be parametrized in 
different ways (see e.g. ||,|57J]). For a narrow width reso- 

11 In perturbation theory a single fermion / of charge Qf and 
color N c f at one-loop contributes 

Aa P = ^ Q2fNcf { i 1 + f ) G{Vf) ~ V! ~ 5/3 } ' (60) 
with 

G{y) = I ^~y ( lo s 1^ - «*) ; < v < 1 

[2y/y- larctan^i=; y > 1 

and j// = Am 2 j / s. A light (£) or heavy (h) fermion yields 



A*? = 







(h) 



The two-loop correction from a lepton is 

Aaf\s) 



~) 



£ + C(3)+ 4 l0g ^f- 1 

+ 3^1og^ + of^ 

s mi \ s A 



(61) 



At least the electron contribution should be taken into account. 
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nance we have 

i? Nw(s) = 9 jM* r W e+e _8(s - Ml) , (64) 

while for a Breit-Wigner resonance 



Rbw{s) 



r R r 



(0) 

R, e+e 



Aa 2 



(V5-M*) a + ^ 



(65) 



We also may consider a field theoretic form of a Breit- 
Wigner resonance obtained by the Dyson summation of 
a massive spin 1 transverse part of the propagator in the 
approximation that the imaginary part of the self-energy 
yields the width by ImiTy(My) = MyFy near resonance. 
Here we have 



Rbw(s) 



9 



r 



(0) 

R, e+e" 



*I% 



Ml 



r B 



(s 



M%)* 



m%*% 



(66) 



Mr and -Tr are the mass and the width of the resonance, 
respectively, and Tr e + e - is the leptonic width as listed 
in the particle data tables. In the formulae above we need 
the undressed leptonic widths M 



r. 



(0) 



R, e+e- 



= r, 



R>e+e - (a/a(M 2 R )) 2 



(67) 



which has been used in this paper: 



F„(a) = 



Ax - A 2 ml 

A 1 + A 2 s -& + f( s ) 
A aa 4 "<; 



s - m 2 + im u r u 



where 



/W = i(^-|) + ia^lo 



2m 7T 



G(s) , (69) 



(i + flr) 



..'ft 



G(s) = 



1 2 j±- 

1 A/ 2 ' M 



(For y/s < m^ + m^ only the real part of G(s) is kept. The 
second term accounts for the p— w-interference. The factor 
G(s) incorporates the effect of the p— ui inelastic channels. 
The parameters are M = 1.2 GeV, r = 0.15 GeV Q and 
n = 0.22, A x = 0.29 GeV 2 , A 2 = -2.3, A 3 = -0.012, 
A4 = 1.84 [pl| .) A slightly modified parametrization (see 
Fig. pL 5| ) has been given more recently with parameters 
fitted to the final CMD-2 data |25J. For other parameteri- 



Analytic formulae for the corresponding real parts the 
reader may find in J57]]. In general one has to take R(s) 
from the data. The imaginary part leads to additional con- 
tributions at order a 2 (s) and in the interference a(s) x 
F 7r (s). In the latter case one has to know the phase of the 
pion form factor, which actually can be determined p3, 
|65|,£36,67|. This issue is beyond the scope of the present 
work. 

In spite of the problems addressed above, what we need is 
the lpi photon VP as a building block for the calculation 
all kind of (e.g., higher order) corrections and this is given 

by 



Fi \s) = [l-Aa(s)}F 7r (s), 



(68) 



which in modulus square agrees with (|13| ) up to contri- 
butions from the imaginary parts. As we have mentioned 
before, in p5[ , corresponding corrections have been ap- 
plied together with the FS correction (|14|) to the "bare" 
cross section referred to as a , v 
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GS fit by CMD-2 2001 
GS fit by Kinoshita 1985 
CMD-2 01 




Fig. 15. The new CMD-2 data on the pion form-factor versus 
Gounaris-Saurai like parametrizations. 



zations see e.g. p3,p3[. In fact the pion form factor can be 
"parametrized" in a more model independent way by ex- 
ploiting analyticity, unitarity and constraints from chiral 
perturbation theory together with information from tttt 
scattering data and by combining iF^js)] 2 data in both 
the space-like and the time-like region |64,p9,p6,67,lL]. 



Pions are not point-like particles. It is therefore not pos- 
sible to calculate the cross sections for pion pair pro- 
duction from first principle using just scalar QED. Usu- 
ally the pion structure is parametrized by the form fac- 
tor F„(s) which contains all non-perturbative QCD effects 
and which is only a function of s. A typical parameteriza- 
tion for F^ is the Gounaris Sakurai parameterization |p9| 



It is the aim to extract F^ from experimental data by un- 
dressing the experimentally observed cross sections from 
radiative corrections. Using the usual procedure of unfol- 
ding the QED corrections leads to a model dependence 
for the results which can be estimated by a comparison 
of different form factor parameterizations. In our radia- 
tive return scenario we can get l-F^s)! 2 directly via (p3[). 
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Note that for this |F ff (M?)| 2 has to be determined at an 
accuracy of 10 % or better. 

Why is the above form factor ansatz a reasonable one 
to parametrize the extended structure of the strongly in- 
teracting bound state pion? First of all it leads to the 
right long range behavior if F^ (0) = 1 which corresponds 
to pure scalar QED. It also allows for a consistent treat- 
ment of radiative corrections under the condition that one 
should not think of a form factor as being related to a pion 
vertex but to the Born amplitude (factorization): 

M (s)[e+e- -> tt+tt-] = M p ° int (s) x F„(a) . (70) 

j^pomt j g ^g B orn amplitude for point-like pions, ob- 
tained from scalar QED. For higher order virtual plus soft 
photon corrections the amplitudes can then be written as 



M v+S (s) = 5 V+S x M p m n \s) x F^s) 

( m i- 

+ terms — ► — 

V 8 



(71) 



The factor 5 V + S is again calculated by scalar QED. Clearly 
this ansatz respects gauge invariance, the renormalization 
procedure of scalar QED can then be applied and the can- 
cellation of infrared divergences is also achieved. 



Note that the "terms^ m\j V stands for non-factorizing 
IFS interference corrections. The above form may be as- 
sumed to hold as an approximation also for the hard pho- 
ton FS corrections. Without further investigations we can- 
not say what is the systematic error we make by utilizing 
this ansatz, however (see the discussion towards the end 
of Sec. |). 
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